Cracklike Dynamics at the Onset of Fractional Sliding 
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We propose an elasto- plastic inspired friction model which incorporates interfacial stiffness. 
Steady state sliding friction is characterized by a generic nonmonotonic behavior, including both 
velocity weakening and strengthening branches. In ID and upon the application of sideway loading, 
we demonstrate the existence of transient cracklike fronts whose velocity is independent of sound 
speed, which we propose to be analogous to the recently discovered slow interfacial rupture fronts. 
Most importantly, the properties of these transient inhomogeneously loaded fronts are determined 
by steady state front solutions at the minimum of the sliding friction law, implying the existence 
of a new velocity scale and a "forbidden gap" of rupture velocities. We highlight the role played by 
interfacial stiffness and supplement our analysis with 2D scaling arguments. 
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The frictional strength and stability of spatially ex- 
tended interfaces is important for a wide range of natural 
and man-made systems [H-@]- Yet, several fundamental 
aspects of it are not well understood. One of these as- 
pects is the onset of frictional sliding in which the (ini- 
tially static) interface separating two bodies in frictional 
contact fails under shear forces, giving rise to relative 
shear motion. 

This transition was recently observed to be mediated 
by the propagation of interfacial fronts ■ By directly 
tracking the real contact area between the two bodies, it 
was demonstrated that in addition to fast sub-Rayleigh 
and super-shear cracklike modes [il 0, [l(|, there exist 
also slow cracklike modes which travel at well-defined 
velocities significantly smaller than the Rayleigh wave 
speed, but nevertheless play a significant role in interface 
weakening. Furthermore, these experiments highlight the 
key role of inhomogcneity to frictional stability and mode 
selection |{|, and the relevance of plastic deformation of 
interlocking asperities to frictional strength These 
important and unexplained observations, especially the 
nature and properties of the slow cracklike fronts, is our 
main focus here. 

In this Letter we develop a dry friction model based 
on the dynamics of microcontacts at frictional inter- 
faces. We demonstrate the existence of a frictional 
instability prior to the onset of sliding, which excites 
cracklike fronts whose velocity is independent of sound 
speed. Most importantly, we show that the proper- 
ties of these slow fronts propagating under transient in- 
homogeneous conditions is determined by steady state 
front solutions at the minimum of the sliding friction 
law, where velocity-weakenin g b ehavior crosses over to 
velocity-strengthening one (3. Our work is strongly 
influenced by recent important numerical investigations 
of a microsco pic ID model [l2| and macroscopic 2D an- 
alytic results 



031 14 



location of the interface and x is the position along it, 

a xy (x,y=0,t) = a(x,t) + r]d t u x {x,y = Q,t) . (1) 

Here a is associated with interfacial contact dynamics, r\ 
is a viscous-friction coefficient and dtu x is the slip rate. 
We choose the velocity-strengthening part r]d t u x to have 
the simplest possible form, which does not affect the gen- 
erality of our results, as long as such a branch exists. In 
addition, we set u y {x, y — 0, t) = 0. 

The next step is to write down dynamic equations for 
<7 in terms of A, the ratio between the real contact area 
and the nominal one, in the form 
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Our starting point is the following expression for the 
total frictional stress a xy (x, y = 0, t), where y = is the 



where A, a xy and u x are functions of x and t at y = 0. 
Eq. ([3]) is analogous to an elasto-plastic decomposition 
e = e tot — e pl (e is a strain measure), corresponding to 
elastic, total and plastic strain rates, respectively. The 
^-function expresses the fact that plasticity is an in- 
trinsically threshold phenomena, where a c is the shear 
strength or yield stress of the microcontacts. The macro- 
scopic stress <j xy is enhanced at the microcontacts level 
by a factor A~ l ~^>\. /iq is the interfacial elastic stiff- 
ness, which is usually not included in friction models (but 
note [HI), h is the effective height of the interface and 
T\ is a basic timescale of irreversible interfacial processes. 
T± 2 = T2 2 + (dtu x /D) 2 , i.e. n is an incoherent sum of an 
intrinsic timescale ti and a kinetic time scale D/d t u x , 
where D is a typical sliding distance [1] . 

The dynamics of A, Eq. consists of two contri- 

butions; the first tends to increase A to a limiting value 
Aq <C 1 on a timescale tq. In general Ao depends on 
the normal stress and may "age" logarithmically on long 
timcscales, but here it is constant. The second contribu- 
tion accounts for the reduction of A due to irreversible 
processes (plastic deformation and eventually fracture) 
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and is similar to the corresponding term in Eq. (j3]). The 
dimensionless parameter k accounts for the possibility 
that A decays on a timescale somewhat different than a, 
as suggested in 0. Eqs. ([T][3]) constitute a "minimal" 
continuum friction model based on the detachment and 
reattachment dynamics of interfacial contacts. A cl osely 
related discrete model has been recently proposed in [12J . 

Focusing on the limit T2^> D /\d t u x \ and defining t- 
t/ro, x=x/D, u = u/D, A = A/A Q , a xy - 
a=a/a c A , Eqs. j21[3]) become 



: & xv /cr c A and 



d- t A = (1 — A) —6 (\a xy \/A - l) KA\d~ t u x \ , (4) 
d~ t d- = aAdiu x - 9 (\a xy \/A- 1J &\dju x \ , (5) 

with a = fj,oD/a c h. Finally, we rewrite Eq. ([1} as 
&xy = d + V^tUx, with f} = r)D / '<7 c tqAq, which completes 
the derivation of our dimensionless friction law. 

Macroscopic frictional phenomena intrinsically involve 
the coupling between bulk elastic deformation and the 
dissipative dynamics of the interface. Therefore, we 
should couple our friction law to a linear-elastic bulk of 
height H. The bulk is described by Lame equation for 
the displacement field u(x,y,t) [l6| . a xy in Eq. (p} is a 
boundary condition at y 
In the limit of small H, a 
condition, but rather a term in a ID equation [1 



(recall that u y (x, 0, t) = 0). 
xy is no longer a boundary 



p %u x = A d xx u x + Od 
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Here p = pHD / '<t c t'qAq and p = pHjo c DA§, where p is 
the mass density and /i is the bulk shear modulus 17 1. 
dd is the dimensionless external driving stress. 

In order to mimic the experimental edge-loading setup 
of [Hi], we study Eqs. for x > with a localized 

driving force of the form ad(x,t) = Kd [vdt—u x (0, £)] 5(x). 
Here Kd and Vd are the rescaled driving spring constant 
and velocity, respectively, and S(-) is a 5-function. 

We first studied our model numerically [l8|. A typical 
solution is presented in Fig. [1] Initially, the response is 
purely elastic, i.e. A = 1 in panel (b) and the driving 
force rises almost linearly in panel (a). At t* (~25 here), 
a cracklike front is initiated at the trailing edge (x = 
0), leaving behind a reduced contact area A < 1. This 
cracklike front is accompanied by a relatively low slip 
rate, see panel (d), but by a significant stress transfer- 
resulting in an inhomogeneous stress state, see panel (c) . 
Finally, this front induces a small, but noticeable, drop 
in the driving force, see the early time deviation from the 
straight dashed line in panel (a). 

Later, at 110), an instability in which a short-lived 
cracklikc front is initiated at the trailing edge and arrests 
at a finite distance (smaller than the system's length L) 
occurs. The front's velocity is much higher than the one 
initiated at t* (compare the slopes in panel (b) and recall 
that the smaller the slope, the higher the velocity) and 



is accompanied by a sharp increase in the slip rate, see 
panel (d). Finally, the transient front transfers a highly 
concentrated stress distribution to the interior of the ma- 
terial, see panel (c), and the jump in slip rate at the edge 
causes an abrupt drop in the driving force, see panel (a). 

After the transient front arrests, the contact area re- 
covers to its background level (in general this will not be 
the case since normal stress variations, which are absent 
here, will carry memory of past deformation), see panel 
(b) , but a highly concentrated stress distribution remains 
inside the material, as a signature of the irreversible slip 
that took place, see panel (c). Another front, similar to 
the one initiated at t*, is initiated at the trailing edge 
soon after (t ~ 130). It eventually triggers another in- 
stability that results in a transient front that propagates 
deeper into the material (not shown). This complex pre- 
cursory activity repeats itself until a front reaches the 
leading edge (x = L), whereupon macroscopic sliding oc- 
curs (not shown), in agreement with experimental obser- 
vations 0- We focus on the first precursory event. 
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FIG. 1: (Color online) A solution of Eqs. (f4I6|) with p = 0.01, 
p = 300, k = 8, a = 1.05, ?7 = 0.1, K d = 2.5 and v d = 0.3. (a) The 
driving force <Xd(t). (b) The contact area A(x,t). (c) a(x,t). 
(d) The slip velocity v(x,t) — dfU x (x,t). For simplicity, the 
tilde's are omitted from the axis labels. 

In order to develop a theoretical understanding of the 
cracklikc dynamics that precedes the onset of sliding, we 
first focus on the elastic response at t < t* . In this regime, 
the contact area remains intact, A= 1, and the 9- function 
in Eqs. (|4][5]) vanishes. Then 



u x (x,i) 



K d v d t 



— exp — x I 9(x) for t < t* 



_(7) 

where the inertial term, which is negligible in this regime, 
was omitted. This solution is valid until a/ A first 
equals unity (note that f]diu x is negligible in this regime, 
thus we can replace a xy with a). This happens when 
&(0,t*) = au x (p,t*) = l, i.e. t* = (K d + y/tfZ)/(aK d v d ). 
The exponential spatial variation in Eq. (JT]) is the ID 
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Green's function (whose analog at higher dimensions is 
a power-law distribution) and hence it characterizes any 
purely elastic region, even at times t>t*. 

We proceed to discuss the cracklike front that initi- 
ates at t*. To understand its nature, we adopt a fracture 
mechanics perspective which tells us that our edge- 
loading system is intrinsically stable. That is, the stress 
level near the tip of a crack is a decreasing function of the 
length of the crack under constant loading at the trailing 
edge. That means that the tip region, which has to satisfy 
some fracture criterion (in our case o xy j A~ 1), cannot 
propagate unless the external loading is increased. This 
implies that the dimensionless velocity c of such a crack 
must satisfy c oc Vd- Our system is actually the frictional 
analog of the famous Obreimoff 's tensile fracture exper- 
iment [l6l | . The prediction c oc Vd is fully supported by 
the numerical solutions. The fact that fracture mechan- 
ics predicts a stable front is important since it immedi- 
ately suggests that the instability taking place at t" has 
a frictional origin, absent in tensile fracture. 

The front that initiates at t* plays an important role 
in triggering the frictional instability at t". In particular, 
it is responsible for the transfer of stress from the trail- 
ing edge to the interior of the material and hence for the 
buildup of an inhomogeneous stress distribution. Such 
inhomogencous stress distributions were shown to play 
an important role in selecting various interfacial cracklikc 
modes in Q . The frictional instability at t~ is accompa- 
nied by a transient cracklike front whose dimensionless 
velocity c>c. Understanding the nature and properties 
of this front is a major goal of this Letter. 

We do not consider here the onset (nucleation) of in- 
stability, which we observed to involve fast, Raylcigh-likc, 
contact area disturbance initiated at the trailing edge, 
but rather focus on the emerging short-lived cracklikc 
fronts. In order to develop a theoretical understanding 
of these fronts we should discuss some basic properties 
of our friction law. Consider the spatially homogeneous 
fixed-points of these equations as a function of the slip 
rate v= d^u x . At "0 = there exist only elastic solutions 
with < a xy < 1. For v > 0, we find the following sliding 
solutions 
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FIG. 2: (Color online) Left: The steady friction law a xy {v) 
(parameters as in Fig. [lj (solid line). The horizontal lines 
correspond to homogeneous loading. Right: The transient 
front velocity c vs. the minimal steady state velocity c™ ~ 
y fia/fjK. The tilde's are omitted from all quantities. 



It is instructive to consider also situations in which 
the homogeneous loading ad, rather than the slip rate v, 
is controlled. In this case, for dd < o c d (see the dashed- 
dotted line in Fig. [2|, we have only one solution (with 
v = 0), which is a stable solution. At dd = a d (see the 
dashed line in Fig. [2]) another solution emerges at v m 
(marked by a filled square in Fig. This allows for the 
propagation of a steady state "phase transition" front 
in which a sliding domain invades a jammed one. For 
&d > o\ we nave two finite v solutions, one with < v < 
v m in the velocity-weakening (unstable) regime and one 
with v > v m in the velocity-strengthening (stable) regime, 
which also allow for steady state front solutions. 

Can one relate the steady state homogeneously loaded 
fronts to the transient fronts observed under inhomoge- 
neous loading? Our basic idea is that the transient fronts 
are short-lived excitations of fronts corresponding to the 
minimum of the sliding friction law. To test this idea, 
denote the minimal velocity of steady state fronts by c™ , 
set dt = -cf s d x in Eqs. (gEJ) and note that Eq. ® im- 
plies v m ~ y/a/fjn and a d ~ f]V m , for k^$>1. Hence, the 
scaling version of Eqs. ([5]) is c™/£ a ~ v m and of Eq. ([5]): 



UV r , 



-fjv m = 



flVr, 



(9) 



1 + KV 



f) v for v > 



(8) 



Since a xy > A for sliding solutions (the ^-function equals 
unity), internal consistency demands that a>l. An im- 
portant feature of the sliding friction law in Eq. ((HJ is its 
nonmonotonic behavior stressed above, where a xy (v) at- 
tains a minimum at v m , with a finite value d' d = a xy (v m ). 
The complete steady state sliding friction law is shown in 
Fig. [2] (left). The parameters a and k, control the rate 
of velocity- weakening for < v < v m , while fj controls the 
rate of velocity-strengthening for v^v m . Their relative 
values determine the crossover at v = v m . We suggest 
that this nonmonotonic behavior is generic [i, 11|. 



Here t s is spatial scale of variation of a, a c d r ~^ f)v m and a ~ 
1 (near the front edge, where the ^-function first equals 
unity) were used and inertia was neglected. Solving for 
and l„, we obtain (with H — H/D) 



\JfLajr\K (10) 



and l a ~ V HI. I = \ijH and l K = fin/ Hr)a > I are two 
rescalcd friction-related lengthscales, and s = fj,/Hfj. 

Eq. (|10p contains an elastic bulk property fj and the 
three friction parameters a, f), k. Adopting our idea that 
c~c™, it provides an analytic prediction for the transient 
front velocity c. We measured c directly from numerical 
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solutions. In Fig. [2] (right) we show c vs. the minimal 
steady state velocity c™ ~ fia/fju. Different symbols 
denote variations of each parameter, shown in the leg- 
end, while the other three were held fixed. The result 
is striking: the data points collapse on a single linear- 
curve with a slope of order unity. This result provides 
strong evidence in favor of our basic idea; the proper- 
ties of the transient inhomogeneously loaded fronts are 
indeed determined by the steady state front solutions at 
the minimum of the sliding friction law. 

A crucial point to note is that Eq. (fTU)) suggests that c 
is independent of the sound speed, c s = yfpjp. For exam- 
ple, in Fig. [TJ c~20 while c s ~170, i.e. c and the sound 
speed are separated by an order of magnitude. Instead 
of p a combination of friction parameters determine c. In 
particular, we highlight the fact that c depends explicitly 
on the interfacial stiffness (1q ~ a UM ■ R is important to 



note that the spectrum of front solutions starts at a fi- 
nite value c™ and hence there exists a "forbidden gap" of 
velocities below it, in contrast to ordinary tensile cracks 
(or standard "phase transition fronts"). This is consis- 
tent with the experiments in 0, 0] , where the slow front 
is always characterized by a well-defined finite velocity. 
Therefore, we propose that the transient fronts observed 
here are analogous to the slow fronts discovered recently 
0, and whose velocity should be described by c™ . 

What about the sub-Rayleigh/super-shear fronts? We 
suspect that this is a matter of the mechanical condi- 
tions at nucleation. When the stored energy levels are 
low, as in the examples considered above, c™ may be the 
most relevant "attractor" for the dynamics. On the other 
hand, when the stored energy levels at nucleation are 
higher, a continuous spectrum of solutions with c > c™ 
may be excited. This seems perfectly consistent with the 
findings of [tj [l9| , where mode selection is shown to be 
controlled by the conditions at nucleation. 

While we have not yet numerically studied our model 
in 2D and hence at the moment cannot comment on the 
validity of our main idea in 2D, we can show that 
varies smoothly from ID to 2D as H is increased and 
hence it is qualitatively independent of dimensionality. 
To see this we note that the ID limit is valid for H-^i£ K ,£. 
Consider then the 2D regime £<^.H. The distinguishing 
feature of dimensions higher than 1 is the existence of 
(cracklike) power-law singularity at £-€i \x\<g.H. One can 
show that in the singular region the boundary conditions 
of the present model reduce to those of 13, |l4j, who 



found that C~ |i| ( - 1+e)/2 , where tan(?re/2) ~c™/s. 



/ln(ff/j. 



If we assume H^>£ K ,£, then our result is identical to 
Eq. (11) in which here reads e~\n(H/£ K 

If we assume £ <S H <g; i K , we can still use the rela- 
tion c™ s j£ s ~ v m , but now with £ s replaced by £ and 
v m by v m (H / £Y 1 ~ e ^ 2 , where the latter describes the 
increase of the slip rate on a scale £ (it decays back 
to v m on a scale H>f). The resulting expression is 



c™/s~J H /£„{£/ 'Hy/ 2 . It matches the ID result of Eq. 

([TU|) , when H ~ £, and the 2D result stated above, when 
H r~j £ K . We reiterate that the most important implication 
of this 2D scaling analysis is that the presence of cracklikc 
singularity does not change the qualitative properties of 
the spectrum of steady state front velocities, which still 
starts at a finite value. 

We believe our results derive from robust properties 
of dry friction and hence are general (and not model- 
specific). In future work we plan to test our ideas in 
2D numerical simulations, to include the effect of normal 
stress variations and aging on the contact area, to ad- 
dress the full spectrum of cracklike fronts and to make 
quantitative comparison to experimental data. 
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